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Abstract. We prove that for any pair of i.i.d. random vectors X, Y in R n 
and any real-valued continuous negative definite function ip : R" — » E 
the inequality 

holds. In particular, for a £ (0, 2] and the Euclidean norm || • |[a one has 

e\\x-y\\% <e\\x + y\\%. 

The latter inequality is due to A. Buja et al. [1] where it is used for some 
applications in multivariate statistics. We show a surprising connection 
with bifractional Brownian motion and provide some related counter- 
examples. 

Mathematics Subject Classification (2010). Primary 60Ef5; Secondary 
60G22, 60E10. 

Keywords. Bifractional Brownian motion, moment inequalities, Bern- 
stein functions, negative definite functions. 



1. Introduction 

Let X, Y be i.i.d. random variables with finite expectations. Then one has 

E\X -Y\ <Epf + Y|. (1.1) 

The inequality (|1.1[) appeared recently in an analytic context (properties of 
integrable functions) [8] . Since (jl.ljl is a nice fact in itself and since it seems 
not to be well known in the probabilistic community, it is desirable to search 
for adequate proofs and to explore possible extensions of it. For instance, for 
which values of a do we have 

M\X — Y\ a E\X + Y\ a ? (1.2) 

As before, we assume that X and Y are i.i.d. and E|AT| Q < oo. 
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Proving (jl.ll) is a non-trivial exercise for a probability course. If X, Y 
are real-valued, one way to see this inequality is to use the identity 

/>oo 

E\X + Y\- E\X - Y\ = 2 / [P(X > r) - F(X < -r)} 2 dr. 

Jo 

For (|1.2p we are, however, not aware of a similar elementary approach. On 
the other hand, A. Buja et al. prove in [I] even a multivariate version of () 1 .2[) : 
for any pair of i.i.d. random vectors X, Y in M. n , any a 6 (0, 2] and for a class 
of norms || ■ || on W l including the Euclidean norm || • H2 the estimate 

E\\X ~Y\\ a ^E\\X + Y\\ a (1.3) 

holds true. The elegance of this inequality is obvious; at the same time we 
stress that it arises from statistical applications. In any case it merits to be 
better known in the probabilistic community! 

In Section [3 we give an extension of () 1 . 3[) by replacing the norm with 
an arbitrary negative definite function. Moreover, we show how this fact 
extends to an arbitrary number of i.i.d. random vectors. In Sections [3J and [¥] 
we establish a surprising connection to some recent advances in the theory of 
random processes related to bifractional Brownian motion. A counterexample 
to (|1.2[) with a G (2, 00) is given in Section [5j 



2. Main result 

Consider the class of continuous real-valued negative definite functions, i.e. 
characteristic exponents of symmetric Levy processes. The notion of negative 
definite function goes back to Schoenberg; good sources are the books [3] 
and [10] . Recall that a continuous real- valued negative definite function is 
uniquely given by its Levy-Khintchine representation 

H0=a+l(Q£,0+ I (l-cos&u))i/(du), £eR n , (2.1) 

Z JR"\{0} 

where a ^ is a constant, Q G R nx ™ is a symmetric positive semidcfinitc 
matrix and v is the Levy measure, i.e. a measure on R™ \ {0} satisfying the 
integrability condition 

/ min{||u||2, \}v{du) < 00. (2.2) 
Jm"\{o} 

Without loss of generality, we will always assume that a — 0, i.e. ip(0) = 0. 
For our discussion it is worth noticing that (£, 77) i-> y — rj) is always a 
metric. A deep theorem of Schoenberg states that a metric space (W 1 , d) can 
be isometrically embedded into an (in general infinite-dimensional) Hilbert 
space H if, and only if, d(£, n) is of the form d^(£, 77) = \/ip(£ — n), cf. [TT] . 
[U p. 187] as well as [7] for a discussion of metric measure spaces related to 
the metric d^,. 
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An important subclass of continuous negative definite functions are the 
spherically symmetric negative definite functions. These are of the form 



Recall that a Bernstein function is a function / : R + — > R + which 
admits the following Levy-Khintchine representation 



here a, b are constants and \x is a measure on (0, oo) satisfying the 
integrability condition J Q min{i, 1} fi(dt) < oo. In probability theory Bern- 
stein functions arise as the characteristic exponents of the Laplace transform 
of subordinators, i.e. increasing one-dimensional Levy processes. Bernstein 
functions, many examples and their connections to various fields of mathe- 
matics are discussed in the monograph [TU|. It is easy to see that Bernstein 
functions are infinitely many times differentiable, increasing, concave; more- 
over, they grow at most linearly. Typical examples are A H> log(l -I- A) and 
A i-> fp(\) := A^ for < (3 ^ 1. Note that the composition / o tp of a Bern- 
stein function / with a continuous real- valued negative definite function ip 
is again a continuous real-valued negative definite function. At the level of 
stochastic processes this corresponds to Bochner's subordination of the Levy 
process with characteristic exponent ip by the subordinator with the Laplace 
exponent /. 

Using the Bernstein functions fp with /3 = a/2 and < a ^ 2 we obtain 



as examples for real-valued continuous negative definite functions. Note that 
the functions defined by (|2.3p are characteristic exponents of subordinate 
Brownian motions. 

We prove the following result extending (|1.3I) . 

Theorem 2.1. Let tp be a real-valued continuous negative definite function on 
M™. For any pair of i.i.d. random vectors X,Y in l n it is true that 



£ ^ /(IIC 1 1 2) where / is a Bernstein function 



(2.3) 




^ IICII2 =/a/2(ll?l| 2 ), 0<a^2, 



Eip(X-Y) ^EiP{X + Y). 



(2.4) 



Proof. Without loss of generality we may assume that a = and Q = - in 
both cases the inequality (|2.4p is elementary. 
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Using the Levy-Khintchine representation of ip we get 
E^{X + Y)=e[ (1 - cos(X + Y, u)) v{du) 

JM n \{0} 

= E / (l -Reexp(i(X + Y,u))) v{du) 

JllL"\{Q} 

(l - Re Eexp(i(X + Y. u})) v{du) 

>\{o} 



>\{0} 



1 - Re[Eexp(i(A, u))f) v{du). 



A similar calculation yields 

EiP(X-Y)=e[ (1 - cos(X - Y, u)) v{du) 

■/R™U0} 



= E / (1 - Reexp(i(X - F,tt))) zy(du) 

JR"\{0} 

= / (l-ReEexp(i(X-l»)Wdu) 
1 - |Eexp(i(X,M))| 2 ) z/(dtt). 



•\{0} 

Using the elementary estimate Re(z 2 ) ^ \z 2 \ = \z\ 2 we obtain (|2.4j) . □ 

Remark 2.2. Let X\, . . . X2 m be i.i.d. random variables in 1" and Sj = ±1 

(non-random, or even random but independent of the X\, . . . , X2 m ) constants 
satisfying J2j=i £ j = 0- Then 

(2m \ / 2m \ 

J^e.xA ^E^ rpxA . (2.5) 

This follows if we use Theorem l2.1l for X — Ylj=i £ t^j an d Y = Ylj=i £ jXj- 



Using the distance function d^(£,,i]) :— -0(£ — rj) related to a real- 
valued continuous negative definite function ip we get the following counter- 
part of (Tj~3|) . 



Corollary 2.3. Let tjj : R™ — > R be a real-valued continuous negative definite 
function, d^^^rf) — y/ — ri) the associated metric and < a ^ 2. For 
any pair of i.i.d. random vectors X, Y in R™ it is true that 

Ed$(X-Y) <E<^(X + Y). (2.6) 

Remark 2.4. Assume that ip : R™ — > R is a continuous function such that 
ip(0) = and = VK - 0- If (|2.4j) holds for this -0 ancl any random 

variable X (and an independent copy Y of X), then one can show that the 
kernel K^(£, rj) := ip(£, + rj) —^{^ — rf) is positive definite. We wonder whether 
this already entails that ip is a continuous negative definite function. 
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3. A relation to random processes 



We will show now that the inequality (|2.4[) has an interesting relation to 
Gaussian processes. Let ip : M™ — > R be a real-valued continuous negative 
definite function defined on M. n . 

Lemma 3.1. The kernel K^(£, rj) = tp(£ + 77) — ip(£ — 77) is positive definite. 
Proof. By the Levy-Khintchine formula (|2.1|) we get 



J^(£, 77) = 2(QC, v) + / (cob(« - »?, u» - cos((e + r?, u») i/(du). 

JE"\{0} 

Using the elementary trigonometric identity 

cos(£ — 7], u) — cos(£ + 77, u) = 2 sin(£, u) sin(?7, u), 

we see that 

K*(£, 77) = 2<Q£, r?) + 2 / sin(C, u) smfo u) i/(du). 

JR"\{0} 

Now let S be a finite set and (A^, £ G 5) be complex numbers. Then 

= 2 V X 6 Xr,(Q^,v) + 2 / J] A e sin(£, u) A,, sin (77, u) i/(du) 

e,7,es •'""XM \tnes J 



\{o} 



2J A 5 sin(£, u) 



v{du) 



>0, 



which means that K^(-, •) is positive definite. □ 

Remark 3.2. A special case of Lemma 13.11 for powers of ^,-norms is proved 
in @]. 

Probabilistic proof of Theorem \2.1\ Since 77) is positive definite, there 

is a centred Gaussian process i e R") whose covariance function is 

For given i.i.d. random vectors X, Y G R™ set 

•= J G|P(dO, 
where P stands for the common distribution of X and Y\ Then 
0<Var(Z*)= /" / K^[C,T])P(d^)P(dT]) 



Mt + ri-W-r,)) P(d£)P(dr,) 
EiJj(X + Y) -E^(X -Y), 
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and we obtain again Eip(X - Y) < Eip(X + Y). □ 

4. Relation to bifractional Brownian motion 

In some most important cases it is possible to identify the Gaussian process 
(Gj, £ £ R n ) of Section [3] with bifractional Brownian motion (bBm). The 
latter process was introduced by Houdre and Villa in [5] as a centred Gaussian 
process B H ' K = (B t ' , t £ R n ) with covariance function 

R"> K {t,s) := E (Bf' K B?> K ) = 2~ K {{\\t\\l H + M^f - II* - s\\ 2 2 " K ) , 

where s,t £ M™. For n = 1, K = 1 we get the usual fractional Brownian 
motion B H with Hurst index H . Originally, the process was defined for the 
parameters H £ (0, 1] and K £ (0, 1]. Bardina and Es-Sebaiy [T] recently 
proved that B H ' K exists for all (H, K) G £>, where 

V :={H,K :Q< H < 1,0 < K ^ 2,H K ^ 1}. 

(The possibility of such an extension was already indicated in the earlier work 
by Lei and Nualart [5] who established an integral representation relating 
B H ' K with fractional Brownian motion B HK ). 
For ip(£) := |£| Q , < a 2, and 

G;f :=2«/ 2 sgn(£)4r> (ei, 
it is trivial to see that 

E(G^'G^) =sgn(e77)2 Q E(4^4r) = ^ + ? 'l"-^-^ a= ^^^- 

Therefore, we are led to a probabilistic interpretation of the inequality (|1.2p 
through Bi' a . 

Remark 4.1. In higher dimensions bi-fractional Brownian motion does not 
show up in the context of our inequalities (nor do we rely on bBm with 
H 7^ h); therefore it becomes natural to search for the extensions of bBm 
based upon general negative definite functions. This will be done elsewhere. 

5. A counterexample 

The inequality (|1.2[) trivially extends to the case a — oo in the following 
sense. Let 

M = sup{r : ¥(X < r) < 1} = ess sup X; 
m = sup{r : P(X < r) = 0} = ess inf X. 

Then 

\\X - YWn - M - m ^ 2max{|M|, \m\} = \\X + Y^. 

Without further assumptions the inequality (|1.2[) will, in general, not 
hold, for 2 < a < oo. To see this, fix a £ (2, oo) and c > 0. For any M ^ c 
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set q := c/M and p := 1 — q. Let Xm,Ym be i.i.d. random variables such 
that 

P(x M = i) = HYm = i)=p; 

F(X M = —M) = F(Y M = -M) = q. 

If M ^ 1, then 

E\X M -Y M \ a -E\X M +Y M \ a 

= 2pq [(M + l) a - (M - l) a ] - 2 a M a q 2 - 2 a p 2 
^ ApqaM"' 1 - 2 a M a q 2 - 2 a p 2 
= M a ~ 2 (4pac - 2 a c 2 ) - 2 a p 2 . 

Hence, whenever c < 2 2 ~ a a and M is large enough, 

E\X M - Y M \ a - E\X M + Y M \ a > 0, 

and (HI fails. 

Remark 5.1. Further counterexamples are presented in 0]. 
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